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A STUDY OF SATELLITE REMOTE MEASUREMENTS OF STRATOSPHERIC 
AEROSOL AND GAS PARAMETERS USING THE SOLAR 
OCCULTATION TECHNIQUE 

By 

William P. Chu^ 

SUMMARY 

This report contains results of research Investigations performed 
during the period from September 1974 to December 1975 under grant 
NSG 1108. Satellite remote measurements of stratospheric aerosol 
and ozone profiles using the solar occultation approach have been 
analyzed based on the known properties of the atmospheric param- 
eters. Atmospheric refraction effects which are important in this 
experiment have been determined analytically. Different inversion 
techniques have been developed based on inverting simulated 
irradiance data. It is shown that aerosol and ozone profiles in 
the stratosphere can be inverted from the irradiance measurements 
to an accuracy of better than 10 percent. 


1 Research Assistant Professor in Physics, School of Sciences and 
Health Professions, Old Dominion University, Norfolk, VA 23508. 


INTRODUCTION 


The purpose of the work performed under this grant was to analyze 
the methods for determining stratospheric aerosol and ozone profiles 
from satellite measurements of atmospheric extinctions during solar 
occultation. Solar occultation experiments have been carried out 
successfully by Pepin (ref. 1) using a high altitude balloon. With 
an orbiting satellite as the observation platform, this technique 
can be used for remotely monitoring stratospheric constituents with 
global coverage. 

In this report we describe the geometry of the experiment and 
derive the basic equation relating the solar radiance to the measured 
irradlance. Atmospheric refraction effects which play an important 
role in determining the geometry of the experiment are analyzed. 

The inversion problems are examined and different inversion 
schemes have been developed. The sensitivity of the different 
inversion techniques to experimental errors has been analyzed. 

DESCRIPTION OF THE EXPERIMENT- 

In the solar occultation experiment, multi-wavelength measure- 
ments of the attenuation of solar radiation along a tangent path 
through the atmosphere are made on board the spacecraft. The basic 
geometry of the experiment is illustrated in figure 1. As the space- 
craft emerges from the Earth's shadow, the radiometer aboard the 
spacecraft will point to the Sun, measuring solar radiation at several 
preselected wavelength bands. As the spacecraft continues in its 
orbit, the line of sight from the spacecraft to the Sun will scan the 
Earth's atmosphere at different tangent altitudes, resulting in a 
measurement of attenuated solar radiation at different atmospheric 
heights. Only two vertical profiles of the atmosphere can be obtained 
for each orbit corresponding to spacecraft sunset and sunrise events. 

Assuming the Earth's atmosphere is divided into several concen- 
tric, homogeneous layers of equal thickness, then the path length 
along the line of sight at a given tangent height is much greater 
for the tangent layer than for any other layers. This implies that 
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the major attenuation will take place within this particular tangent 
layer, and can result in a very high vertical resolution measurement 
of atmospheric constituents. 

EXTINCTION MODEL OF THE STRATOSPHERE 

In the wavelength range from 0.3 to 1.0 micron, stratospheric 
attenuations are predominantly caused by aerosol, ozone, and air 
molecules (Rayleigh component) with minor contributions from NO 2 , 
water vapor and O 2 molecules. The peak of the aerosol distribution 
generally occurs at an altitude of 18 km while the peak of the ozone 
profile occurs at 23 km. Models of the stratospheric extinction 
versus wavelength profiles can be constructed from different vertical 
distributions of aerosol and ozone. Figure 2 shows an extinction 
coefficient model as a function of wavelength at the aerosol peak 
altitude of 18 km. Figure 3 shows the extinction model at the ozone 
peak altitude of 23 km. The aerosol vertical distribution used in 
generating these plots was based on Fenn's (ref. 2) moderate volcanic 
model, while the ozone profile was based on Elterman's model (ref. 3). 
The aerosol extinction as a function of wavelength can be calculated 
from the particle size distribution, total density, and index of 
refraction. Figure 4 shows the same extinction model at an altitude 
of 30 km. NO 2 extinction which was computed from recent measurements 
(ref. 4) shows some contribution to the total extinction at high 
altitude. 


BASIC EQUATIONS 

The equations relating the observed irradiance to the atmospheric 

extinction profiles can be derived by considering the geometry as 

shown in figure 1. For a radiometer situated on a satellite with 

orbital velocity vector v, a single ray originating from the Sun 

with intensity I , after propagating through the atmosphere with 

s 

tangent height h,_ will reach the radioriieter making an angle 6 
with respect to the vector v. The intensity reaching the radiometer 
1(0) will be given by 
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1 ( 6 ) 




( 1 ) 


where T^(h,p) Is the transmission function of the atmosphere at 
wavelength \ and is given by the Lambert-Beer Law as 


Tj^(h^) 



6j^(h) dPj^(h) 


( 2 ) 


3^(h) is the total extinction profile of the atmosphere at wavelength 
X including contributions from aerosol, Rayleigh component, and other 
gases. is the ray path length at wavelength X. 

Assuming the atmosphere is spherically stratified and spherically 
symmetric, then each sunray making a view angle of 6 with respect to 
the vector v at the radiometer will define a tangent height h^(6). 
The solar radiance distribution Nj^(0) at the entrance window to the 
radiometer will be given by the following equation 

N^(e,t) - Sj^(6,t) T^[h^(6)] (3) 

Sj^(0,t) is the refracted but unattenuated solar radiance profile 
at time t which is determined by the solar limb darkening profile at 
wavelength X. The effects due to atmospheric refraction will be 
discussed in more detail in the next section. Both and 

are in units of Watts/cm^-£2. The irradiance received by the radiom- 
eter will then by given by 

*^X^®i^ " /^X^®'®i^ Sj^(0,t^) Tj^[h^(0)] dfl (4) 

where Fj^(0,0^) is the radiometer’s field of view (FOV) function 
centered at view angle 0^, and Hj^(0) is the irradiance received 
by the radiometer at time t^ when the FOV is pointing at view angle 

“i- 

Equation (4) is the desired result relating the observed irradi- 
ance to the source function and the transmission function of the 
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atmosphere The transmission function is independent of 

the spacecraft geometry and is defined solely by the extinction pro- 
file 8(h) of the atmosphere through equation (2). 

ATMOSPHERIC REFRACTION 

Atmospheric refraction is an important factor to be considered 
in the solar occultation experiment. The effects can be dramatic, 
especially in the lower stratospheric layers where bending of sunrays 
is pronounced. In the following section we will describe a ray trace 
method to calculate atmospheric refraction effects. 

1. Ray Trace Analysis 

Consider a spherically symmetric atmosphere, the equation govern- 
ing a ray path is (ref. 5) 

n(r) r sin * constant (5) 

where n(r) is the refractive index at position r with respect to 

the Earth's center, 0 is the angle between the position vector and 

the ray tangent vector at r. Referring to figure 5, we assume the 

spacecraft is situated at point P at a distance of r from the 

o o 

Earth's center. Defining Cartesian x-y axes at the Earth's center, 
we start a ray propagating at an angle 9^ from the x-axis. We 
have from equation (5) 


n(r ) r sin A « C 
' o' o ^o o 


where r_ * (x^ + y^)^/^, and <]>^ is related to 0 as 
o o o o o 


(P^ a COS“^ 




( 6 ) 


(7) 


where i is the unit vector along the x-axis. We now trace a small 
segment of ray with length As and arrive at point Pj at a distance 
rj from the Earth's center where 
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ri - j (x^ + As cos e^)2 + (y^ + As sin 0^)2 j 


1/2 


( 8 ) 


The ray direction at can now be computed from equation (5) 


to be 


(^x * sin 


-1 ( ^2 \ 

\n(ri)ri / 


( 9 ) 


and the new angle 0x is given by 


01 - cos-1 ( ^ ^ - ♦ 


( 10 ) 


^e can continue this process until the ray passes through the 
atmosphere and determine the final angle 0^. Moreover, the tangent 
height of *:his ray can be determined from the r value when ^ ■ ir/2. 

The above procedure can be done for a large number of rays with 
different Initial ray angles. We then generate a table of three 
para.meters; the initial ray angle 0^, the final ray angle 0^, 
and the corresponding ray tangent height h^. The relationship 
between these three parameters will completely describe the refrac- 
tion effects. 

2. Ray Trace Results 

The ray trace method has been applied to atmospheres with dif- 
ferent temperature and pressure distributions. The refractive index 
n(r) of air depends on its density as (ref. 6) 

n(r) - 1 + p(r) (11) 

where p(r) is the density of air determined by the temperature and 
pressure distributions, and is a constant which differs for 

different wavelengths. 

Figure 6 shows results from the ray trace calculation for a 
satellite at orbital altitude of 600 km viewing the Sun. An orbital 
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period of 90 minutes was used to compute the 
Sun shape for a sunset event. The atmosphere used was the standard 
1962 mid-latitude atmospheric profile. The top and bottom of the 
solar disk as observed on the image plane of the radiometer are 
shown as the satellite heading toward the Earth's shadow. It is 
obvious that the Sun images are continuously being demagnlfled. 


Figure 7 shows the computed Sun Images as observed on the space- 
craft for a sunset event. The vertical dimension of the Sun images 
can be demagnlfled by a factor of four as the apparent height of the 
Sun is at 10 km altitude. Computations of the refracted but unattenu- 
ated solar radiance profiles S^(6,t) can now be performed based on 
the refraction calculation as derived above. 


INVERSION TECHNIQUES 

In the solar occultation experiment, the radiometer aboard the 
spacecraft will measure solar irradlance as the line of sight of the 
Instrument scans across the vertical tangent height of the atmosphere. 
Equation (4) will reduce to a set of m equations describing the 
irradlance measured at m different view angles 


Hj^(Bi) - Jfj^( 6,0^) Sj^(0,t) exp [- 2 /* 6j^(h)dP(h)j 

h^(6) 


dn 


( 12 ) 


1*1, 2, . . ,m 


Given a set of measurements equation (12) is nonlinear 

in the unknown B^(h). The techniques used for the solution of equa- 
tion (12) will depend on the functional form of the instrument FOV 
function Fj^(0,0^). We will discuss this in the following two 
special cases of interest. 

1. Case 1: ^ Radiometer with Small FOV 


If the vertical resolution of the radiometer's FOV is small com- 
pared to the vertical structure of the atmosphere which is of interest 
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we can assume the transmission function x^[h,^(6}] to be constant 
for each irradiance measurement. Defining an effective optical 
depth we have from equation (2) 


gf - sT $j^(h) dP(h) (13) 

\(6i) 

where is related to the measurement by 

- in| f F^(6,e^) Sj^(e,t^)dn/H^(e^)j (U) 


Equation (13) is a Fredholm integral equation of the first kind. 
The integral can be approximated with n discrete atmospheric layers 
with equal thickness, and the equations can be replaced by the follow- 
ing linear system of equations 


g^ - E Pij 1 "* 2,...m (15) 


where 

layer, 

height 

matrix 

is the total extinction coefficient 
and is the ray path length in the 
at the ith layer. Equation (15) can be 
form 

in the Jth atmospheric 
jth layer with tangent 
abbreviated into a 


g - P6 

(16) 

Equation (16) can be solved directly for 
when n m or in the least-square sense when 

6 either exactly 
n < m. 


$ - P"ig n ■ m 

(17) 


e - (p'^?)-ip'^g n < ra 

(18) 


T 

where is the inverse of P, and P is the transpose of P* 
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The direct inversion schemes are generally unstable due to the 
presence of noise term associated with the measurements g^. In 
or'^'^r to suppress unwanted oscillations in the solutions, the 
sm. thing method developed by Twomey (ref. 7) and Phillips (ref. 8) 
can be used. The smoothing method Imposes certain constraint on 
the solutions, and the most widely used constraint has been the 
minimization of the second-difference expression 
*»■ resulting solutions are given by the following form 

6 - (P*^? + yH)“» p'^g (19) 

where H is the second-difference constraint matrix and y is the 
smoothing parameter whose magnitude is dependent on the noise level 
in the measurements. 

2. Case 2: Radiometer with Large FOV 

This case generally does not arise unless the radiometer is 
measuring the irradlance from the whole Sun. The inversion of equa- 
tion (12) in this case generally begins by linearizing the integral 
equations in some ways. One form of linearization uses an expansion 
of the solution around some initial guessed profile 

e^(h) - S®(h) + AB3^(h) (20) 

Assuming that 


2^ ABj^(h)dP(h) « 1 (21) 

hip 

we can reduce the transmission function to the following form, 

Tj^lh^(6)l - T°[h^(0)l [l - 2J AB;^(h)dP(h)} 

h^ 


( 22 ) 


where 


T° Ihj(e)) - exp j-2y 8® (h)dP(h)j (23) 


Substituting equation (22) into equation (12), and assuming 
■ 1, we obtain 

Ch,p(e)] 

(24) 

X A8;^(h)dP(h)j dn 


where 

hJOj) -rs^(e,t^)Tj [hj(e)]dO 


(25) 


The double integral in equation (24) can be approximated by the 
following system of equations 


AH. 


n n 

L E 

J-l k-J 


SlJ •’jk 


(26) 


where AH^^ » Hj^(0^) - H^(6^), and S^j is the solar irradiance 
contribution in the ith measurement from the jth ati.'ospheric layer. 
Defining ®ij ©Quation (26) then becomes 


■ k?i 


(27) 
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or in abbreviated matrix form 


LH - KAB (28) 

Equation (28) is identical in form to equation (16) so that the 
previously discussed methods of solution can be applied here. However, 
the initial guessed profile this case could be quite poor. 

A method of solution in which the initial guessed profile is con- 
tinuously updated can be done with the iterative scheme. One for»:\ 
of the iterative scheme which is based on Twomey's work (refs. 7, 9) 
is 

3^+^ - + k'^(Kk'^ + Y)"l - H^J (29) 

m 

0 

whe.* e 6 is the tth order inverted B profile, is the measured 

0 ® 
iri.vdiance, H is the tth order computed Irradlance, and y is the 

ratio of noise variance to signal variance. The iteration will stop 
when the difference between the computed irradlance and the measure- 
ments are within the inst^rument noise level. Another form of the 
iterative scheme which was developed by Conrath and Revah (ref. 9) 
and Smith (ref. 10) is 

» 3^ + G-l k'^(H„ - H*") (30) 

m 

where G is a diagonal matrix with elements * 

comparison of the two iterative schemes in satellite remote tempera- 
ture probing haL been discussed by Conrath and Revah (ref. 9). 

SENSITIVITY ANALYSIS 

The sensitivity of the inversion techniques to various experi- 
mental errors has been analyzed based on the inversion of different 
computer simulated irradlance profiles. These analyses have been 
performed in order to establish the accuracy of the inversion 
approaches and the accuracy requirement for different instrumenta- 
tion design parameters. 
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The first step in generating the expected irradiance profiles 
is to compute the atmospheric transmission function several 

prese- 3 Cted wavelengths. We have chosen wavelengths at 0.38, 0.45, 

0.6, and 1.0 ym to coincide with the spectral channel locations for 
the Stratospheric Aerosol and Gas Experiment (SAGE). The 0.6 micron 
channel is the ozone channel while the other three are primarily 
aerosol channels. Figure 8 shows the calculated transmission versus 
altitude profiles for the four wavelengths based on Elterman’s aerosol 
and ozone models (ref. 3). These calculations have included atmos- 
pheric refraction effects based on the standard 1962 mid-latitude 
atmospheric temperature and pressure profiles. 

To incorporate the spacecraft geometry in generating the irradi- 
ance profiles, we consider a radiometer with field of view of 30 arc 
seconds, situated on a satellite at orbital altitude of 600 km. The 
radiometer is assumed to point at the solar disk at 10 arc minutes 
from the vertical top edge during the complete sunrise or sunset event. 

The inversion of these irradiance profiles was then performed 
by introducing experimental errors of various origins and magnitudes 
to the irradiance profiles. There are two types of errors which will 
directly affect the inversion results. The first type of error is 
noise associated with the desired signals. This random noise has 
both optical and electrical origins. The optical contributions are 
dominated by irradiance detected from outside the instrumental spectral 
bandwidth and outside the instantaneous field of view. The electrical 
contributions are the usual quantization error in the A to D converter 
and .noise in the electronics. The second type of error is associated 
with the pointing system. Since the radiometer's pointing can be con- 
trolled by a servo mechanism in maintaining a fixed location on the 
solar image, a slight fluctuation in the servo system can produce 
uncertainty in the instantaneous viewing direction and result in 
tangent-height error. These tangent-height errors when considering 
different amounts of attenuation at different tangent altitudes, plus 
the solar limb darkening profiles, will generate irradiance errors. 

The two types of errors mentioned above in different magnitudes 
were added to the expected irradiance profiles and inversions were 
performed to determine their effects. The atmosphere was divided 
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into homogeneous layers of 1-km thickness with inversion performed 
for each layer. Figures 9 through 12 show the direct inversion 
results based on equation (17) for the four channels simulated 
irradiance profiles. An i/istrumental noise level of 1 percent 
of reading and pointing errors of 3 arc seconds (one sigma) have 
been included. The results indicated that aerosol extinction pro- 
files from the three aerosol channels can be inverted with an accuracy 
of about 10 percent from 10 to 30 km altitude. Similary, ozone 
extinction profiles can be extracted to about 10 percent accuracy 
after the aerosol contribution has been subtracted. 

The large oscillations in the direct inversion results can be 
suppressed using the smoothing scheme as described by equation (19). 
Inversions were performed for the same irradiance profiles using the 
smoothing scheme with a second-difference constraint on the solutions. 
The smoothing parameter y used in the inversions is inversely pro- 
portional to the logarithm of the effective transmission T. The 
reason for choosing a nonlinear variable y can be seen by referring 
to equation (14). The measurements g^|^ are related to the effective 
tranmission T^ as 

g^ = In (T^) (31) 

Irradiance errors in the measurements can be included in the 
quantities g^ as 

g^ = In (H^ + AH^)/H^ (32) 

where AH^ are the errors associated with the irradiance measurements 
H. , and H is the unattenuated irradiance. For AH. to be small 
compared to d we can expand equation (32) and obtain 

g^ = ln(H^/H^) + AH^/H^ (33) 


The first term on the right side of equation (33) is the signal 
term, while the second term is the noise term. In choosing the y 


to be inversely proportional to the signal to noise ratio, we then 
arrive at the particular form of y as discussed above. 

Figures 13 through 16 show the inversion results with the 
smoothing scheme. Considerable Improvement in the solutions has been 
obtained with the fine structure oscillations being suppressed. 

Aerosol extinction profiles can be reproduced to better than 10 per- 
cent accuracy, while ozone extinction profiles can be reproduced to 
better than 5 percent. 

Other models for the vertical distributions of aerosol also have 
been used in simulating the four channels irradiance profiles for 
inversion. Figures 17 through 19 illustrate the inversion results 
from the smoothing scheme for the three aerosol channels based on 
Fenn’s background aerosol model (ref. 2), while figures 20 through 
22 illustrate inversion results for a moderate volcanic aerosol dis- 
tribution (ref. 2). The noise levels in these figures are relative 
to the full scale irradiance levels. The accuracies of the inverted 
results for different aerosol distributions are generally similar. 

CONCLUSIONS 

In this study we have examined the retrieval of the vertical 
distributions of stratospheric aerosol and ozone from multi-wavelength 
measurements of transmitted radiation in the solar occultation exper- 
iment. Inversions of simulated irradiance profiles computed from 
different aerosol and ozone vertical distributions have been performed 
at 0.38, 0.45, 0.6, and 1.0 ym wavelengths. The inversion results 
indicated that aerosol extinction profiles can be retrieved to an 
accuracy of better than 10 percent over most of the stratosphere, 
while ozone vertical distribution can be retrieved to better than 
5 percent accuracy. 

The accuracy of the inversion techniques as discussed in this 
study is not dependent on the method of measurements. Irradiance 
data can be obtained either with the radiometer pointing to a fixed 
location on the solar image during the occultation event or with the 
radiometer scanning across the solar image. Provided that the 
irradiance errors associated with each measurement are within the 
same level, accuracy of the retreived profiles will be similar. 
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This study has been concerned mainly with the retrieval of 
multi-wavelength extinction profiles for the stratospheric consti- 
tuents. Since different wavelength extinction coefficients of the 
stratospheric aerosol can be measured at each altitude, these 
wavelength-dependent extinction coefficients can be used to retrieve 
some physical properties of the stratospheric aerosol such as size 
distribution, number density, and index of refraction. A future 
study to evaluate this problem will be very important. 




I 



REFERENCES 

1. T.J. Pepin, Ph.D. thesis, University of Minnesota (1970). 

2. R.W. Fenn, unpublished works. 

3. L. Elterman, AFCRL-68-0153 (1968). 

4. J. Fontanella, A. Girard, L. Gramont, and N. Loulsnard, 

Applied Optics, 14, 825 (1975). 

5. M. Born, and E. Wolf, Principles of Optics, Pergamon Press, 

Oxford, 1964, p. 123. 

6. P.N. Tverskoi, Physics of the Atmosphere, Israel Program for 

Scientific Translations, p. 457 (1965). 

7. S. Twomey, J. Assn. Computing Machinery, 10, 99 (1963). 

8. D.L. Phillips, J. Assn. Computing Machinery, 9, 97 (1962). 

9. B.J. Conrath, and I. Revah, NASA TM X-62, 150 (1972). 

10. W.L. Smith, Applied Optics, 9, 1993 (1970). 


16 



Figure 1. Solar occult at ion geometry 





Figure 2. Extinction coefficients model as a function of wavelength 
at an altitude of 18 km. 
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Figure 7. Sun shape images as observed on the spacecraft for different apparent 
altitudes from the refraction calculations. 
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Figure 10. Direct inversion results from simulated irradiance profile at 0.45 micron 
wavelength. 
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Figure 17. Inversion results with smoothing at 1.0 micron for a background aerosol 
model. 
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Figure 18, Inversion results with smoothing at 0.45 micron for a background 
aerosol model. 
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